Contact Angle Effects on Pore and Corner Arc Menisci in Polygonal Capillary Tubes Studied with the Pseudopotential Multiphase Lattice Boltzmann Model

In porous media, pore geometry and wettability are determinant factors for capillary flow in drainage or imbibition. Pores are often considered as cylindrical tubes in analytical or computational studies. Such simplification prevents the capture of phenomena occurring in pore corners. Considering the corners of pores is crucial to realistically study capillary flow and to accurately estimate liquid distribution, degree of saturation and dynamic liquid behavior in pores and in porous media. In this study, capillary flow in polygonal tubes is studied with the Shan-Chen pseudopotential multiphase lattice Boltzmann model (LBM). The LB model is first validated through a contact angle test and a capillary intrusion test. Then capillary rise in square and triangular tubes is simulated and the pore meniscus height is investigated as a function of contact angle θ. Also, the occurrence of fluid in the tube corners, referred to as corner arc menisci, is studied in terms of curvature versus degree of saturation. In polygonal capillary tubes, the number of sides leads to a critical contact angle θc which is known as a key parameter for the existence of the two configurations. LBM succeeds in simulating the formation of a pore meniscus at θ > θc or the occurrence of corner arc menisci at θ < θc. The curvature of corner arc menisci is known to decrease with increasing saturation and decreasing contact angle as described by the Mayer and Stoewe-Princen (MS-P) theory. We obtain simulation results that are in good qualitative and quantitative agreement with the analytical solutions in terms of height of pore meniscus versus contact angle and curvature of corner arc menisci versus saturation degree. LBM is a suitable and promising tool for a better understanding of the complicated phenomena of multiphase flow in porous media.